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Abstract. The gravitational potential and the gravitational rotation field generated by an thin-disk mass distribu- 
tion with exponential density are considered in the case when the force between any two mass elements is not the 
usual Newtonian one, but some general central force. We use an approximation such that in the Newtonian case the 
gravitational field generated by the disk reduces to the familiar expression that results from applying the Gauss' law. 
In this approximation, we invert the usual integral relations in such a way that the elemental interaction (between two 
point-like masses) is obtained as a function of the overall gravitational field (the one generated by the distribution). 
Thus, we have a direct way for testing whether it is possible or not to find a correction to the Newtonian law of gravity 
that can explain the observed dynamics in spiral galaxies without dark matter. 
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.,— i 1. Introduction. 

X' ■— i 

$H In a previous paper (Rodrigo-Blanco 1996 : hereafter, paper I ) we found the exact solution to the problem of, given 
. the gravitational field generated by an spherical mass distribution with exponential density, finding out which was 
the elemental gravitational force, that is, the force between any two point-like masses, that could have generated that 
field. 

In this paper we follow the same steps as in paper I , but the mass distribution that we consider is a thin disk with 
exponential density, which is a more realistic description of a spiral galaxy. 

In this case, the exact solution cannot be found and the problem is solved in an approximation that, in the 
Newtonian case, is equivalent to use the Gauss' law for calculating the gravitational field. We show that this is a very 
good approximation which improves even more when the g c s(f), that represents the deviation from the Newtonian 
force, is a growing function of r, which is a welcome behaviour for explaining the observed data. We also see that this 
approximation is clearly better than using the results obtained in paper I for the case of spherical symmetry. 

In Section we give the general definitions that are used later in Section || for the case of a thin disk with exponential 
density. In Section [I| we show the mathematical basis underlying the results presented in Section ||[ Finally we offer 
some conclusions. In Appendix [A| we study the validity of the approximation. In Appendix |b| we list some of the 
mathematical identities used in Section [|. 



2 



C.Rodrigo-Blanco: Non-Newtonian gravity directly from observations: thin-disk mass distribution. 



2. General definitions 

Let us assume that the gravitational potential generated by a point-like mass does not correspond to the usual 
Newtonian form but can be written in terms of a function g(r) that describes the deviation from the Newtonian law, 
that is, 

Hr) = -^9(r). (1) 
The force per unit mass is, by definition, the gradient of the potential, 

, Gam , . r , . 

F « = -^-9eff(r)-, (2) 

where we have introduced 

g cS (r) = g{r) - rg'(r). (3) 

In this way, to find the total potential or the total force generated by a mass distribution £1 with density p(r), it 
is necessary integrate over the volume spanned by f2 to get: 



F(R) = -Go / / / S r "jj*-;" „ (r > (5 ) 



for the potential experienced by a point mass at a distance R from the centre of and 

s(\R-r\) R-r 
\R~r\ 2 \R-r\ 

for the force. 

In the case that the gravitational potential is only a function of the distance to the centre of the distribution, it is 
convenient to introduce two new functions ^f(R) and ^ e s(R) such that: 

•(*) - - ^f* *(*), (6) 
F(R) = * cS (R) f (7) 

and the rotation velocity of a test particle in a circular orbit bound to the distribution will be: 

vLiR) - ^f^ *e ff (i?) (8) 
where the auxiliary functions ^ c s(R) and ^(R) satisfy the following relationship: 

* eS (R) = *{R) - RV'(R), (9) 

Our goal is to design a procedure where, assuming that F(R) is known (say from observation of the rotation 
velocity) for all values of R, a g c g (r) that generates the given rotation velocity is obtained. In other words, given the 
potential as inferred from observations we want to find which g(r) could have generated it. Actually, we find g{r) and 
g e fi(r) as functions of ^> c s(R) and ^{R) respectively. 

3. Thin-disk mass distribution with exponential density 

The luminosity profile of many spiral galaxies can be well fitted assuming that the luminous matter is placed along a 



thin disk with a density that decreases exponentially with the distance to the centre of the galaxy (Kent 1987) 



p(r) po e- ar . (10) 

Considering a thin-disk distribution and using Eqs. (||), (0) and (|Io|), in Eqs. ^ and (||), the two problems outlined 
above can be recast as two integral equations: 
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i) Given ^(R), denned by 

2tt 



find a function g(r) such that: 



dr I dd 
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r2 — 2Rrcost 



\/R 2 + r 2 - 2Rrcos( 
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or 



R 



(11) 



(ii) Given ^ e s(R), defined by ([?]), find a function g c s(r) such that: 



dr / dO 

n Jo Jo 



2tv 



g cS WR 2 +r 2 - 2Rrcos( 
(R 2 +r 2 - 2Rrcos6)§ 



■(R — r cos 9)r 2 sin( 



R 2 ' 



(12) 



In this case the problem cannot be solved exactly. We use an approximation that we call Gaussian as, in the 
Newtonian case, it is equivalent to use the Gauss' law for calculating the gravitational field. 

In the next section the calculations are described in detail. The solution to the two problems outlined above, in 
the Gaussian approximation can be summarized as: 
(i) Potential problem: (viz. Eqs. (0) and (11)) 



In this case, the approximate solution to the problem is 

g(x) = *(x) - \*"(x), 
or 

where the function \1/ has the following behaviour at the origin: 
*(0) = 0. 

(ii) Force and velocity problem: (viz. Eqs. ([?]), (||) and ([l2])). 

Here, the approximate solution is given by the following expression: 

g ctf (x) = V cS (x) - + 4-*eff(z) 

a z a z x 
and the behaviour of ^ at the origin is as follows: 

* cff (0) = <ff(0) = 0. 



(13) 
(14) 

(15) 
(16) 



4. Mathematical formalism. 



First we show how to solve the potential problem, that is, how to go from equation (|llj) to equations (13) and (t#. 
Then, the results obtained above are used to tackle the force problem, i.e, to go from equation ( p^ ) to equations (15) 
and (|l|). 

4-1. The potential. 

In order to go from equation Eq. ( pi] ) to Eqs. (|l^) and (jl4|), it is necessary to use several mathematical identities. For 
convenience, these are listed in Appendix [3. 



Using ([Bl]), (|B3|),(|B|) and (|B7|) into (|ll|) *(i?) can be written as: 

2 



r(2fc+i/2) 
r(2fe+i) 



io°° d P S£ ]T J 2k+i/2(pR) J °° dre ar \ZrJ 2k+ i/ 2 (pr) 



(17) 



We apply our approximation here. We only consider the first term in the series of Bessel functions, that is, we 
drop all the terms in the series but the one with k = 0. It can be seen, and we prove it in an appendix, that in the 
Newtonian limit (i.e, when g(r) = g c s(f) = 1) this approximation corresponds to apply the Gauss' law to the mass 
distribution, i.e, to say that the gravitational force at a distance R is proportional to the total mass contained in the 
sphere of radius R. This is the reason why we have called this method Gaussian approximation. As it is quite a good 
approximation in the Newtonian case, and we also will see that it is exact for some forms of g(r) and a very good 
approximation for the most interesting forms of g(r), we restrict ourselves to consider only the first term in the series. 

Thus, we will write: 



MR) = - 



dp^J l/2 ( P R) 
o P JO 



dre ar ^/r.J l i 2 (pr) 



(18) 
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Now, using the functional form of J\/2 and applying (B2) to invert the Fourier transformation, and after some 
straightforward calculations, we get: 



*o(i2) 



„ct-R 



dxg(x)e 



dxg{x)e a0L 



-aR 



R />oc 

dxg(x)e ax — / dxg(x)e~ 
o Jo 



(19) 



(20) 



At this point, it is useful to introduce an auxiliary function ip(x) that makes the integrals exact: 

g(x) = ifj{x) - —^i/j"(x). 

or 

Eq (|2^) can be inserted into Eq. (jl^) and, upon integration by parts, we get: 

* (i?) = V(i?) - e~ ai V(0), (21) 

where ip is a solution to the ordinary differential Eq. (|2^) that satisfies the conditions of being an analytic function at 
x = 0, and 



lim ip(x) exp(— ax) = lim ip'( x ) ex P( — ax ) = 0. 



(22) 



These conditions are easily satisfied in all the astrophysical systems. The analiticity is satisfied in the Newtonian 
limit, which is the behaviour that we expect to recover at R ~ and it can be seen that if Eq. ( p2| ) were not 
satisfied, the rotation velocity would grow almost exponentially with the distance, which clearly seems to contradict 
the observations. 

It is straightforward to see that, provided tp(R) is a solution to Eq. (|20|), then ^q(R) is also a solution to the same 
equation. Moreover, the term proportional to ip(0) in Eq. (|l|) implies that \&o is zer o at the origin. Taking all that 
into account, we finally obtain: 



g(x) = tf (a:) 



1 



2 *"(*), 



and 

*(0) = 0. 

where we have omitted the subscript everywhere as we will in what follows. 



(23) 



(24) 



4-2. The force and the rotation velocity. 

Once Vl'(-R) is known, ^> s(R) can be calculated using Eq. (||). Equivalently, once g(r) is known, g e s(r) can be obtained 
through Eq. (||). Using these two equations together with Eq. (p3|), and after some straightforward calculations, a direct 
relation between g e g and \I/ e g- can be obtained: 



g ctf {x) = y cff (x) - \^ s {x) + -^-*eff(z)- 



(25) 



From Eq. (^) and the behaviour of 4" at the origin, Eq. (p4j), it is easy to see that, at the origin, ^ e ff will satisfy: 
* cff (0) = <f' cS (0) = 0. (26) 



5. Summary and conclusions. 

We have found the solution to the problem of inverting the integral relation between the elemental law of gravity and 
the overall gravitational field generated by a thin-disk mass distribution. We have done it in an approximation that 
we have called Gaussian as it is equivalent to use the Gauss' law for calculating the gravitational field generated by 
the distribution. Although this is not exact, we show in an appendix that it is a very good approximation, and it gets 
much better when g c g (r) grows with r, which is the expected behaviour if the observed gravitational behaviour must 
be explained without the need of dark matter. 

In summary, we now have a method for inferring g e s(f) given 4' c ff(i?). It can be said that we have a way to travel 
from the world of macroscopic interactions to the world of microscopic interactions. It can be schematized as follows: 
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Given the observed V ro t{R) for a given galaxy, use Eq. (g) to obtain ^> e ff(R), fit it by a mathematical function and 
then use Eq. (jl^) to get the g e ff(r) that describes the elemental gravitational force (through Eq. (ph) that can explain 
the observed rotation curve. 

We are now ready for applying the differential expression that we have found to the observed rotation curves of 
spiral galaxies, and find the g e ff required for explaining those rotation curves. Doing that, we will see whether it is 
possible or not to find a universal law of gravity that can explain all the rotation curves needin g only the observed 



luminous matter. This will be done in a separate publication (Rodrigo-Blanco & Perez-Mercadcr 1996 ) 
A. The validity of the approximation. 



In the step from Eq. ( pL 7[ ) to Eq. (18) we made the approximation of considering only the first term of an infinite 
series. The terms dropped depend also on g(r), which is not known, so it is not trivial to evaluate the goodness of the 
approximation. 

In this appendix, we choose some functions g c s(r) and compare the exact rotation velocity with the approximate 
one. By exact we mean the rotation velocity generated by the considered g e g(r) through Eqs. ( |l2| ) and (||) . Nevertheless, 
for most functional forms of g c s(r) an analytical solution cannot be found and it is necessary to perform numerical 
integration for finding that exact rotation velocity. Whenever that happens we have chosen to perform the integrals 
assuming that the disk has some non-zero thickness. Actually this is a more realistic model for a spiral galaxy, the 
astrophysical system to which our approximate method is applied. By approximate rotation velocity we mean the one 
such that the corresponding \l/ c g(i?) satisfies Eqs. ( f25[) and (|26|). 

As a test of consistency we first consider the Newtonian case, i.e., g(r) — g c s(r) = 1. When we use Eqs. (^5|) and 
( ^6|) we obtain: 

* off (i?) = l-(l + ai?)e- Qfl , (Al) 
and thus, the rotation velocity is 

KldR) = - (1 + aR)e-« R ] = (A2) 

where M(R) is the disk mass inside the sphere of radius R. Of course, this is not the exact result, but it is what we find 
if we apply the Gauss' law as an approximation for evaluating the gravitational field. That is why the approximation 
is called Gaussian. 

Next, our approximation must be checked for other different forms of g e g. For doing that we choose a parametric 
family of <7 e ff(?")'s given by 

g cS ^(r)={£j\ (A3) 

where /i parametrizes how fast g c s grows. 

Concerning the exact solution we have found the analytical one (for the thin-disk case) for [i = 1 and ji = 3: 

VL^(R) = - (1 + oK)e-«] = * (A4) 

^3(«> = ^(f) a (l + _^). m 

For the other values of \i the integrals are performed numerically for a disk with a small thickness h — a^ 1 /6 
Concerning the approximation, it can be seen that, for /i = 3 the solution is (|A5| ), that is, the approximation is 
exact. 

In Figure || we have plotted the exact rotation velocities compared with those corresponding to our approximation 
for some values of \i. It can be seen that the results are very similar in every case. For each case, the solutions are 
normalized dividing by a convenient constant Vq defined as: 

G M tot a . . 

Vo = —, — r— • (A6) 
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Fig. 1. Exact rotation curves (dots) and Gaussian approximation (solid line) when g e ft(r) = (^j for some values of fi (/i — is 
the Newtonian case). In every case, for the sake of clarity, the velocities are normalized by dividing by an appropriate constant 
Vo = G °*^'° . In the inset graphs (dashed lines) we have plotted, as a function of r for each case, the percentage of error made 
when the Gaussian approximation is used instead of the numerical integrals. 




Fig. 2. Mean square error in the velocity when the Gaussian approximation is used instead of the exact solution for the disk, 
as a function of /i, and using g e a(r) = The exact meaning of o"d is explained in the text. 

Although it can be seen in Figure |l| that the approximation is quite good in every case, and is better when g e g(r) 
is a growing function of r, it is interesting to have a more quantitative way of describing the difference in the rotation 
velocity obtained in the exact case and in the Gaussian approximation as a function of fi. In order to do that, we 
define a quantity er^ as follows: 

°°M = nZ vim ■ (A7) 

where the sub-scripts D and G stand for the disk (exact) and the Gaussian approximation, respectively. We sum over 
r.i, which are the points where the integrals are calculated. The total number of points for each value of /x is N = 100. 
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So defined, <r|)(/i) is a measure of the mean square error that we make in the rotation velocity if the approximation 
is used instead of the numerical integrals, for each value of p. 

In figure (0), we plot the value of a 2 D versus p. Once again, it can be seen that the larger p is, the smaller the 
difference. We also plot the values of aj-, obtained in paper I when the equations for spherical symmetry were used 
as an approximation to the disk problem. It can be seen that the Gaussian approximation is always better than the 
spherical one. 

B. Appendix: Mathematical identities 

1. Fourier sine transform 

g(\fr 2 + R 2 - 2rRcos8) = - / g s (p) sm(p^/r 2 + R 2 - 2rR cos 9)dp (Bl) 

n Jo 

9s(p) = / g(t)smptdt (B2) 



2. Addition theorem for Bessel functions (see Gradshteyn 1980) 



= Hl» f> + k) J " +k(mp) J - +k(mr) C£(cos 9) (B3) 



k=0 



where: 



lo = \Jr 2 + p 2 — 2rp cos 9 

p<r (B4) 
= Gegenbauer Polynomials. 



Using Eq. (B3), together with: 



sin(mz) = ^^J 1/2 (mz) 



we get: 

sin p\J r 2 + R 2 — 2rR cos ( 



Vr 2 + R 2 -2rRcos8 2*jRr f^L 

k— 



(B5) 



(2k + l)J k+ Upr)J k+i (pR)P k (cos9) (B6) 



3. Legendre Polynomials. 

f 2 [n^0)] 2 ;if fc = 2n + l 
d9P k (cos9) = < (B7) 

[o ; if k = 2n 
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